

df <- data.frame(
  Roll_1 = sample(1:6, size = 1000, replace = TRUE),
  Roll_2 = sample(1:6, size = 1000, replace = TRUE)
)

t <- df %>% 
  mutate(Sum = Roll_1 + Roll_2) %>%
  filter(Roll_1 > 3) %>%
  group_by(Sum) %>%
  summarise(Count = n()) %>%
  mutate(Probability = Count / sum(Count))


135 / 553
